X-ray variability of AGNs in the soft and the hard X-ray bands 
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ABSTRACT 

We investigate the X-ray variability characteristics of hard X-ray selected AGNs (based on 
Swift/BAT data) in the soft X-ray band using the RXTE/ASM data. The uncertainties involved 
in the individual dwell measurements of ASM are critically examined and a method is developed 
to combine a large number of dwells with appropriate error propagation to derive long duration 
flux measurements (greater than 10 days). We also provide a general prescription to estimate the 
errors in variability derived from rms values from unequally spaced data. Though the derived 
variability for individual sources are not of very high significance, we find that, in general, the 
soft X-ray variability is higher than those in hard X-rays and the variability strengths decrease 
with energy for the diverse classes of AGN. We also examine the strength of variability as a 
function of the break time scale in the power density spectrum (derived from the estimated mass 
and bolometric luminosity of the sources) and find that the data are consistent with the idea of 
higher variability at time scales longer than the break time scale. 

Subject headings: black hole physics — galaxies: active — galaxies: nuclei — galaxies: Seyfert — X-rays: 
galaxies 
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1. Introduction 

Time variable X-ray emission is a probe to the 
inner regions of the Active Galactic Nuclei (AGN) 
and is considered to be one of their defining char- 
acteristics. X-ray variability studies provide in- 
sight into the geometry and the physical condi- 
tions in the nuclear regions on account of the fact 
that this emission is thought to be emitted from 
regions close to the supermassive black hole. The 
X-ray variability studies on time-scales of months 
to years in AGNs have established the similarity 
of the physical processe s across stella r mass to su- 
permassive black holes (|Uttlev fc McHardv . 200i 



McHardv et al. 20061 ). Since the characteristic 



time-scale is proportional to the mass of the black 
hole, Mbhi continuous monitoring for a long dura- 
tion becomes necessary in the case of AGN based 
on linear scaling with black hole mass from X-ray 
binary systems. 

Early studies on X-ray variability using data 
from EXOSAT in the 0.1 - 10 keV range showed 
that on short time-scales AGN variability ap- 
peared to be red noise dominated. In other 
word s, it was unpredictable and aperiodic in na- 
ture ( McHardv fc Czernv Ill987t ). Corresponding 
power spectral density (PSD, that is, variability 
power as a function of temporal frequency) is best 
fitted by a power-law of slopes -1 to -2 with no cut- 
off seen down to the lowest sampled frequencies. 
However, the shapes of their PSDs were shown 
to be similar to t hose of X-ray bin aries (XRBs) 
in their soft state ([McHardv 1119881 ). Hence anal- 
ogous to the high-frequency breaks seen in the 
PSDs of XRBs and assuming the break time-scale 
varying linearly with Mbh, a break to flatter 
PSD slope was expected over a period of days to 
weeks. Due to the uneven sampling of data the 
results yielding the break-frequency were uncer- 
tain. The EXOSAT data of a ~day-long AGN 
X-ray observations have revealed inverse corre- 
lation between t he amplitude of variability and 
X-ray luminosity (jBarr fc Mushotzkv |[l986[ ). Re- 
sults from the ROSAT All-Sky Survey (RASS) 
on the soft X-ray variability of AGN showed the 
variability strength on timescales of days to be a 
function of steepness of the X-ray spectrum with 
sources with steeper spectra exhibiting stronger 



variability ( Grupe et al. |[2001 ). PSD studies were 
constrained due to large time gaps between sub- 
sequent observations in ROSAT data,. 

With the launch of the RXTE in 1995, a sig- 
nificant improvement was seen in the quality of 
data and a manifold increase in the monitor- 
ing timescale, from months to years, was possi- 
ble. Because of its rapid slewing capability and 
flexible scheduling, an evenly sampled long-term 
monitoring of AGN X-ray variability was car- 
ried out for several sources. Contrary to the re- 
sults from the then existing satellites, AGN light 
curves from RXTE/PCA in the 2-20 keV energy 
band showed that on longer timescales, of about 
a month, sources displayed less dispersion in vari- 
ability amplitudes c ompared to those measured o n 
time scales of 1 day ( Markowitz fc Edelson |[2001 ) . 
The PSDs also revealed a cutoff/ break at long 
timescales correlated with t heir black hol e mass 
('Edelso n fc Nandra I 1 19991 : lUttlev et al~l 120021: 
Marko witz et al. 1 12003). but with a l arge scat- 
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ter in the correl a tion (IDone fc Gierlinski I |200i). 
McHardv etlTl (|2006l ) ascribed this scatter to 
a third variable, accretion rate, and found that 
the break time scale combined with accretion rate 
can predict the mass of the black hole all the 
way from Galactic X-ray binaries to supermas- 
sive AGNs. They postulated that the variability 
originates within the accretion disk and the break 
time scale is associated with the inner edge of the 
disk, which makes an inward movement with an 
increasing accretion rate for a given black hole 
mass. 

Since such a scenario makes a definitive sug- 
gestion of the accretion disk geometry, a study of 
variability with energy should be able to pin down 
the radiation processes in the inner accretion disk 
quite reliably. Recently, the variability studies of 
AGN at energies a bove 15 keV are carried out by 
Beckmann et al. I (|2007t ). The data of the first 9 
months of the Swift /B AT aW-sky survey in the 14 - 
195 keV range of the 44 brightest AGN revealed a 
tendency of unabsorbed or type 1 Seyfert galaxies 
to be less variable than absorbed or type 2 objects. 
Also they found a more solid anti-correlation be- 
tween variability and luminosity, which was pre- 
viously detected in soft X-rays, UV and optical 
bands. 

The RXTE/ ASM all-sky survey data gives in- 
formation on a large number of X-ray sources, 
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including AGNs. The sensitivity (typically 10 
mCrab for the one day average data) is not suffi- 
cient to make a detailed study of AGNs. By tak- 
ing data at larger bin sizes, however, it should be 
possible to get meaningful information on the vari- 
ability at longer time scales. In this paper, we have 
carried out a comparative study of the variability 
of AGN in soft (< 12 keV) and hard (> 12 keV) 
X-ray bands using data from the RXTE/ASM all- 
sky survey and results obtained by Beckmann et 
al. from S'wi/i/S^T all- sky survey. Data selection 
is described in section 2 followed by data analysis 
and results in section 3 and discussion and conclu- 
sions in section 4. 

2. Data Selection 



The All-Sky Monitor (ASM, iLevine et aT 



( 19961 )) onboard Ross i X-ray Timing Explorer 
(i^Xm ls^alkl (|l999f )) consists of three Shadow 
Scanning Cameras (SSCs). Each SSC contains 
a position-sensitive proportional counter (PSPC) 
that views sky through a slit mask. FOV of each 
SSC is 6°x90°, allowing ASM to scan most of the 
sky every 1.5 hours. So apart from locating tran- 
sient objects, it also provides photometric records 
of known sources in three energy bands corre- 
sponding to A (1.5 - 3 keV), B (3 - 5 keV) and 
C (5 - 12 keV) in addition to the total, or sum 
band, intensity in the 1.5 - 12 keV band. The 
MIT database (http://xte.mit.edu/ASM_lc.html), 
from where we have extracted ASM light curves, 
gives data dwell by dwell and one day average data 
points. In dwell by dwell each raw data point rep- 
resents the fitted source flux from one 90 second 
dwell. Whereas in case of 'one day average' each 
data point represents the one-day average of the 
fitted source fluxes from a number (typically 5 - 
10) of individual ASM dwells. 

To have a comprehensive information, we have 
selected all AGN candidates (147) from the ASM 
source list and carried out a systematic analy- 
sis .^_We_Jouiid_3]^AGN_from the sample used 
bv iBeckmann et al. ( 2007 ). We find measurable 
ASM flux for all the sources except two (NGC 
1365 and GRS 1734-292). This sample of com- 
mon sources consists of 3 blazars, 1 radio galaxy, 
8 Seyfert 1, 4 Seyfert 1.5, 2 Seyfert 1.8, 1 Seyfert 
1.9 and 12 Seyfert 2. The BAT sources not in the 
ASM source list include 1 blazar, 3 Seyfert 1, 1 



Seyfert 1.5 and 8 Seyfert 2. 

3. Data Analysis and Results 

3.1. ASM dwell data selection criteria 

We have downloaded all available ASM dwell 



data for 147 AGN from MIT database ( |http://xte 
mit.edu/ASM_lc.html). This covers data from 
MJD 50087 to about 54466, i.e. 1996 January 5 
to 2008 July 15. Dwell data are subjected to the 
usual selection cuts as prescribed by the ASM web- 
site and binned in 20 and 40 days bins. On closer 
inspection of light curves of individual sources, 
several large spikes, some of them with an appar- 
ent periodicity of one year, were noticed. These 
abnormal data were found to have large intrinsic 
errors too. Since the count rates in ASM is based 
on a profile fitting method, the errors depend on 
the other sources in the field of view and the ori- 
entation of the ASM FOV for a given observation. 
As a first cut, we examined the possibility that 
there could be enhanced unaccounted systematic 
errors for data points with large measured errors. 
It is found that the measured errors for each dwell 
data is sharply peaked at 1 counts/s. The distri- 
bution of errors for each dwell data for one of the 
sources (IC4329A) is shown in Fig. [TJ Based on 
these considerations, we have introduced a selec- 
tion criterion of an upper limit on sum count rate 
error of 3 counts/s, corresponding to thrice the 
peak value of the distribution. 

After selecting dwell count rates for sum, A, B 
and C bands based on the above mentioned cuts, 
they are binned into 20 and 40 days starting from 
MJD 50087. For each bin the average rate and er- 
ror are calculated for all the four bands, which are 
weighted mean and estimated error of the mean. 
Data points with too few number of individual 
dwell measurements (less than 20 data points for 
20 day binning and less than 40 data points for 40 
day binning) are also ignored, ensuring that there 
is at least one dwell per day, on an average. For 
example, for IC4329A, Fig. [2] shows the sum band 
rate error after binning as a function of number 
of dwells per bin, for 20 and 40 days binning. In 
both the cases sum band error increases sharply 
for bins with lower number of dwells. After intro- 
ducing these criteria, the individual light curves 
are carefully examined and no systematic abnor- 
malities are found in these light curves. 
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To summarize, the selection criteria for gener- 
ating ASM light curves are: 



IC4329A dwell data 



^ 2000 - 




sum bond roti 



Fig. 1. — Distribution of errors on sum band rates 
from dwells for IC4329A 



1. Reduced chi-square of the fit < 1.5 

2. Number of sources in the field of view < 16 

3. Earth angle > 75° 

4. Exposure time > 30 seconds 

5. Long-axis angle -41.5° < 6 < 46° 

6. Short-axis angle -5° < < 5° 

7. Error on sum count rate < 3 counts/s 

8. Number of dwells for N day bin light curve 

is > N 

The first six selection criteria are the same as 
those used for generating one-day average count 
rates from dwell data as described in the MIT 
database webpage. 



IC4329A : 20 doys binning 




IC4329A : 40 days binning 



Fig. 2. — Sum rate error vs no. of dwells for 
IC4329A with 20 days and 40 days binning. 



3.2. ASM light curves 

The ASM measurements have a systematic 
error of 3% added in quadrature to the sta- 
tistical errors. This error is estimated using 
Crab light curves and could be underestimate 
in some cases ( http: / / xte.mit.edu/ AS M_lc.html| ) . 
Grimm et all (|2002f ) have studied light curves of 
several sources with constant X-ray flux binned 
with different bin durations from 1 to 200 days 
and have estimated systematic error on large time 
scales in the range of 0.01-0.1 cts/s, depending 
on the source flux. According to Fig. 2 of their 
paper, there seems to be a shift of about 10% be- 
tween measured and estimated RMS for sources 
with constant flux. We have investigated this in 
detail to get a better handle on the average sys- 
tematic errors. For a few of the bright sources, we 
have generated the power spectral density spectra 
(PSD) using a method based on the autocorrela- 
tion function (Gilfanov & Arefiev 2005) and have 
assumed that the frequency independent power 
has to be due to the errors in each observations. 
Based on these PSD, rms noise level seems to be 
underestimated by a factor of 1.13. A few steady 
sources too were investigated and the frequency 
independent systematic error is similar to this 
number. Though we cannot completely rule out 
some additional frequency dependent systematic 
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errors, we assume that for large bin sizes (greater 
than 10 days), most of the systematic errors are 
taken into account. This assumption is further 
strengthened by examining individual light curves 
as well as the distribution of the source fluxes (see 
below). After correcting typical binned error of 
0.005 c/s with this factor, this systematic error 
translates into typical systematic error of 0.5 c/s 
for dwell in S band and 0.3 c/s for A, B and C 
bands. 

According to ASM data products page (http:// 
xte.mit.edu/), there is a 1 mCrab positive bias 
in the light curve intensities. This bias is only 
evident when enough data are averaged so that 
the statistical uncertainties are driven to be very 
small. This generally requires binning or analyz- 
ing the light curves on time scales of many days 
or longer. In order to estimate this bias we have 
carried out the following exercise. We have gen- 
erated histogram of average sum band count rates 
for 147 AGN (Fig. [S]). Unlike the usual number 
distribution given in an integral way, this is a dif- 
ferential plot and hence the number should vary 
with count rates as a power-law with an index of - 
2.5 for an isotropic distribution. Hence we have 
fitted this histogram with a function consisting 
of a power-law with slope -2.5 and constant off- 
set represented as the start point of the power- 
law. Count rates were spread in the X-axis using 
a Gaussian distribution with (7=0.0056 counts/s 
according to the above mentioned error estima- 
tion. Fit for sum band data is shown in the Fig. [3] 
The number distribution is shown as data points 
with error bars (error taken is the square root of 
the number in each bin) and the fit is shown as 
a histogram. Offsets obtained with this method 
for various bands are as follows: S band : 0.0895 
c/s, A band : 0.0274 c/s, B band : 0.015 c/s and 
C band : 0.0213 c/s. These values are consistent 
with the 1 mCrab offset quoted in the ASM web- 
site. 

As a further verification of these offsets and 
systematic error, we checked for consistency of 
flux estimates. We compared our results with the 
flux obtained from an uniform sp ectral analysis 



done on a large sample of AGNs (jWinter et al 



20091) and found 22 comm on sources. Us ii ig the 



( 20091) 



spectral parameters from IWinter et al. 
and ASM sum band average count rates we 
estimated 2-10 keV absorbed flux for these 



sources using web calculator (WebPIMMS tool at 
Ihttp: / /heasarc . gsfc . nasa.gov /Tools ) . Fig. |4] shows 
the plot of ten years averaged absorbed ASM flux 
vs observed flux from other experiments, both 
in the energy range of 2-10 keV. It should be 
noted that these observations are not simulta- 
neous. Hence source variability can cause non- 
identical fluxes between the ASM and the other 
experiments. Points are marked by diamonds. 
Solid line corresponds to the case where ASM 
flux is equal to flux measurements from other ex- 
periments. Dotted lines correspond to ASM flux 
equivalent to ±0.1 c/s around solid line. Most 
of the sources barring a few like IC4329A and 
NGC 4388 are close to the fitted fine, within ±0.1 
ASM c/s. As a further check we compared the 
long duration RXTE PCA measurements from the 
uniform analysis done by iMarkowitz fc Edelson 
(|2004f ). We found these measurements for nine 
sources: IC4329A, NGC4151, NGC3783, MCG-6- 
30-15, 3C1 20, 3C390.3, NGC3227' NGC 4051 and 
NGC3516 (|Markowitz fc Edelson I I200I . We esti- 
mated the ASM flux by extracting 40 days binned 
sum band rates near simultaneous with the PCA 
measurements. These points are indicated by plus 
signs with error bars in the plot. Hence overall 
there is a good agreement between flux estimates 
from ASM and other observations. This agree- 
ment also justifles our estimates of count rate 
offsets and systematic errors. We note here that 
further reflnement in the flux comparison is diffi- 
cult to achieve because of the fact that at shorter 
time scales (where pointed observations are avail- 
able) ASM error bars are too large and at longer 
time scales (when the ASM error bars are driven 
to a small value) pointed observations are simply 
not available. 

Applying the above mentioned cuts and correc- 
tions, light curves in all four bands were generated 
for all ASM detected AGNs. Average sum band 
count rates of all these AGN in four bands with 20 
days binning for 30 sources common to ASM and 
BAT database are given in Table [T] There was 
one more common source, GRS 1734-292, which 
was rejected after applying above the mentioned 
cuts. That is, there was not a single stretch of 20 
days with at least 20 dwells satisfying all the cuts 
(possibly because of its proximity to the Galactic 
center). Count rates given in Table [T] are weighted 
mean of binned counts rates and estimated error 
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0.1 0.2 0.3 

Count Rate (1.5-12 keV) 



Fig. 3. — Histogram of ASM Sum band count rates 
for 147 AGN. Data is shown as points with error 
bars and the fit by power law with constant offset 
is shown as a histogram. 
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2,0x10 '4.0x10 '6,0x10 ''8,0x10 '1,0x10 ^1,2x10 '^,4x10 
Flux from other experiments (2-10 keV) 



Fig. 4. — Absorbed ASM flux vs observed flux 
(in units of ergs cm~^ s~^) from various exper- 
iments shown by diamonds. ASM flux is de- 
rived from average of ASM count rate over the 
period of more than ten years. Sohd line cor- 
responds to to the case where ASM flux is the 
same as flux from other measurements. Dot- 
ted lines correspond to ASM flux in the range of 
±0.1 count/s around solid hne. Points indicated 
by plus signs (and large error bars) correspond 
to flux measurements from PCA onboard RXTE 
compared with simultaneous ASM measurements 
for 1C4329A, NGC4151, NGC3783, MCG-6-30- 
15, 3C120, 3C390.3, NGC3227, NGC4051 and 
NGC3516. 



on mean. Out of these, 15 sources have average 
sum band count rates above 0.2 counts/s. Some of 
these sources (11) are common sources in BAT and 
ASM database and are included in Table [TJ The 
remaining sources are lES1959-h650, IGR J18027- 
1455, MCG 6-30-15 and Mrk 501. Light curves for 
these 15 sources with 20 days binning are given in 
Figs. El Eland m 

The literature survey for RXTE-PCA Hght 
curves of AGNs (in this study) accumulated over 
long periods has allowed us to compare them with 
those of long looks of ASM. We could obtain PCA 
light curves for Mrk 501 and Mrk 421 taken over 
a period of 8 and 9 years respectively overlapping 
with the period for ASM light curves. The data set 
for Mrk 421 covers a period 1996-2005 (Figure 1 of 



Emmanoulopoulos fc Wagner 



jioO?)) and 1997- 

lUboe). Peaks at 



2004 for Mrk 501 ()Gliozzi et al. 
similar MJDs are identified on inspecting the ASM 
and PCA light curves for Mrk 421. The peaks were 
found near MJDs - 51000, 51700, 51900, 52600 and 
53100. Similarly, peaks for Mrk 501 were found 
near the time 1997 ( MJD 50600). 

3.3. Variability strength 

We have calculated the strength of variability 
for all these sources in bin sizes of 20 and 40 days. 
The strength of variability and the errors on them 
are calculated as follows. 

Light curve from source consists of N flux mea- 
surements Xi with measurement errors of ct^. In 
addition to these variations, object has intrinsic 
variability or additional source variance ag. It is 
necessary to disentangle these two variances. One 
of the approaches to estimate the intrinsic variabil- 
i ty is to use an e xcess variance axs as an estimator 



(jNandra et al. 111997 : .Vaughan et al. 2003) . It is 



given by 



^2 



where the sample variance S^ is given by 



N 



X? 



(1) 



(2) 



where x is mean rate and af is average variance 
of the measurements. 

For light curves with varying measurement un- 
certainties {ffi 7^ constant), it is necessary to 
use a numerical approach to obtain the best es- 
timate for the parameter of interest (ctq, here) 
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Table 1 

Variability strength and Tereafc for ASM sample 



Source 


ASM rate 


ASM 


ASM variability strength^ 


BAT^ 


BH Mass 


log(L6oi) 


log(TB) 


Type 


name 


(1.2-15 keV) 
counts 


luminosity 




S„ 




variability 


(log Mbh) 


ergs s~^ 


days 






log(L^SM) 
ergs 








strength 


Msun 












1.5-12 keV 


1.5-3 keV 


3-12 keV 


S„ 










3C273 


0.280±0.004 


45.68 


34.7±3.4 


64.7±7.9 


33.6±4.4 


15±5 


7.22 


47.35 


-3.04 


Blazar 


3C454.3 


0.078±0.006 


46.59 








42±12 


9.17 


47.27 


1.13 


Blazar 


Mrk421 


0.894±0.005 


44.57 


83.7±4.2 


79.4±4.1 


88.4±4.5 


142±38 


8.29 


~45.00 


1.51 


Blazar 


IGR J21247+5058 


0.126±0.004 


43.54 


42.3±5.0 




51.4±5.0 


11±6 








Radio galaxy 


3C120 


0.173±0.005 


44.11 


42.7±6.3 


64.1±7.9 


38.9±10.4 


<10 


7.42 


45.34 


-0.65 


Syl 


3C390.3 


0.109±0.004 


44.37 


29.7±5.3 


49.6±9.4 


30.5±7.5 


<6 


8.55 


44.88 


2.17 


Syl 


EXO0556-386 


0.069±0.005 


43.74 




76.8±14.1 




<9 








Syl 


IC4329A 


0.459±0.006 


43.91 


22.7it2.5 


48.0±4.5 


18.1±3.5 


<3 


6.77 


44.78 


-1.47 


Syl 


MR2251-178 


0.134±0.005 


44.57 


51.2±8.4 


96.7±22.5 


42.3±13.5 


<7 








Syl 


NGC3783 


0.222±0.005 


43.16 


28.5±4.9 


67.3±12.8 


24.1±7.4 


<4 


6.94 


44.41 


-0.75 


Syl 


NGC4593 


0.105±0.005 


42.76 


79.7±11.2 




71.3±21.8 


<7 


6.91 


44.09 


-0.50 


Syl 


NGC3227 


0.125±0.005 


42.10 


56.2±7.1 


123.2±17.6 


49.3±11.4 


<14 


7.64 


43.86 


1.26 


Syl. 5 


NGC3516 


0.091±0.004 


42.68 


76.6±7.0 




72.0±9.6 


<7 


7.36 


44.29 


0.25 


Syl. 5 


NGC4051 


0.064±0.005 


41.31 








<9 


6.13 


43.56 


-1.62 


Syl. 5 


NGC4151 


0.434±0.005 


42.51 


53.9±3.2 


136.3±12.5 


53.3±3.2 


27±7 


7.13 


43.73 


0.32 


Syl. 5 


NGC1365 


<0.005 










<17 








Syl. 8 


MCG-05-23-016 


0.263±0.005 


43.11 


32.8±3.6 




32.2±4.6 


6±4 


7.60 


44.21 


2.85 


Syl. 9 


NGC5506 


0.237±0.006 


42.79 


42.9±4.6 




47.7±5.2 


<3 


8.00 


44.53 


1.36 


Syl. 9 


CenA 


0.674±0.006 


42.19 


32.7±1.9 


73.4±7.9 


35.9±2.0 


10±2 


8.38 


43.0 


3.66 


Sy2 


CygnA 


0.331±0.005 


44.85 


25.8±2.4 


36.0±4.1 


27.7±3.0 


<7 


9.40 


46.0 


2.86 


Sy2 


ESO103-G35 


0.120±0.006 


43.16 


oy.oity.u 


Tft S 4- 1 /I 1 


1 4.0It iU. 1 


<8 








Sy2 


Mrk348 


0.054±0.006 


42.91 








12±10 


7.21 


44.27 


-0.04 


Sy2 


NGC1275 


2.541±0.005 


44.73 


4.62±0.41 


5.6±0.8 


4.7±0.5 


<14 


8.51 


45.04 


1.93 


Sy2 


NGC2110 


0.168±0.005 


42.84 


41.7±5.3 


69.7±13.9 


28.4±7.6 


25±7 


8.30 


44.10 


2.41 


Sy2 


NGC2992 


0.122±0.005 


42.69 


101.7±8.1 


128.1±14.4 


118.6±11.1 


45±19 


7.72 


43.92 


1.37 


Sy2 


NGC3081 


0.022±0.005 


41.98 








23±11 








Sy2 


NGC4388 


0.204±0.004 


42.99 


17.6±7.8 


36.0±10.6 


32.6±15.3 


11±4 








Sy2 


NGC4507 


0.062±0.006 


42.77 








<8 








Sy2 


NGC7172 


0.063±0.006 


42.52 








12±9 








Sy2 


NGC7582 


0.056±0.005 


42.02 








23±21 








Sy2 


Mrk3 












±7 


8.65 
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0.254±0.010 
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42.3±5.0 
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70.4±12.1 










Syl 


MCG-6-30-15 


0.208±0.006 


44.65 


34.3±5.7 


69.1±7.8 


40.2±8.7 




6.65 


43.56 


-0.52 


Syl 


Mrk 501 


0.368±0.004 


43.27 


76.7±3.9 


62.1±3.6 


90.4±4.7 




9.21 






Blazar 


Crab 


75.355±0.010 




0.48±0.03 


1.25±0.07 


0.37±0.03 


1.27 








Pulsar 



^ASM variability strengths computed for data spanning duration of 12.5 years, whereas BAT variability strengths correspond to data collected over 9 months. Bin size 
is 20 days for both the cases. 
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Fig. 5.— ASM Sum band light curves (1.5-20 keV) 
for lES 1959+650, 3C 273, Cen A, Cyg A and IC 
4329A with 20 days binning. 



( Almaini et al. ] I2OOOI) . The most widely used 
method for this purpose is the principle of max- 
imum likelihood. The probability density of ob- 
taining data values Xi is given as a product of 
Gaussian functions. Using Bayes' theorem proba- 
bility distribution for ctq is obtained. This is like- 
lihood function for <tq and it can be calculated as- 
suming Bayesian prior distribution for erg and xt. 
By differentiating, maximum li kelihood estimate 
can be obtained for ctq (see 



Beckmann et al. 



(j2007l ) for equations). In the case of identical mea- 
surement errors {ai = constant) this expression 
reduces to excess variance for uniform prior. This 
corresponds to aQ=axs- In the present analysis 
of ASM data, since ai is almost constant, we have 
used this simplified approach. This assumption 
about approximate constancy of aiS is established 
using IC4329A data. Fig. [8] shows the distribu- 
tion of sum rate error after 20 days and 40 days 
binning, after applying cut on no. of dwells per 
bin. Mean value of binned sum rate error and 
RMS is given in the figure. These distributions 
are sufficiently narrow to validate the assumption 
of constancy of diS. 

Variability is given in terms of normalized ex- 
cess variance, i.e. ^nxs~'^xs/^'^ '-'^ the frac- 
tional root mean square (rms) variability ampli- 
tude (Fyar) given by 



IGR J18027-1.|55 



MCG-5-23-16( 



MCG-6-30-15 



Mrk 501 



5x10" S.lxlO'' 5.2>^W 5.3x10" 

Time (in MJD) 



Fig. 6.— ASM Sum band light curves (1.5-20 keV) 
for IGR J18027-1455, MCG-5-23-16, MCG-6-30- 
15, Mrk 412 and Mrk 501 with 20 days binning. 



p — 



x^ 



(3) 



Error on variability or normalized excess variance 
consists of two parts: 1. arising from measurement 
error and 2. arising from intrinsic fluctuations, 
depending on the index of the power spectrum. 
So err{al,xs) is given by 



[err{aljxs)\T 



[err{alixs)]li ' 



[err{<Tj^xs)?i 
(4) 

where subscript T stands for total error, M for 
measurement error and I for intrinsic error. Ac- 



cording to equation (11) from IVaughan et al. 
(I2OO3I) . 



[err{a'^i^xs)fM = 



2 {alf , 4afFi 



N 



N 



x^ 



(5) 



Also approximating intrinsic fluctuations with 
white noise (it should be noted that the light 
curves are red noise and white noise is assumed 
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here fo r simplicity) and following I Vaughan et al, 
(120031 ). 
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Fig. 7.— ASM Sum band light curves (1.5-20 keV) 
for NGC 1275, NGC 3783, NGC 4151, NGC 4388 
and NGC 5506 with 20 days binning. 



IC4.329A : 20 days binning 

Meon = 0-0962 +/- 0.0485 J 




Sum bond rote error 
IC4329A : 40 days binning 



Mean = 0.0666 +/- 0.0324 



0.05 0.10 0.15 0,20 

Sum band rate error 



Fig. 8. — Distribution of sum rate errors with 20 
days and 40 days binning for IC4329A after cut 
on number of dwells. 



Combining equations 5 and 6 we get 



Substituting for Fyar from equation 3 we get 

2 



[err{a%xs)]T 



Hence 



N [ x^ 



(6) 



(7) 



(8) 



(9) 



Following equation (B2) of I Vaughan et al. I(|2003l ). 

1 



err{Fvar)T 



IF 



-pM<^nxs) (10) 



This is valid when si^ficrj^xs) small and we get 
the following expression 



err{F^arh 



1 



/2Nx^F^, 



(11) 



When err{a'ffxs) l^-rge, the error estimate will 
only be approximate. The strength of variability is 
calculated from the expression Sy — lOO%.aQ/x 
i.e. Sv = lOO%.Fyar. 

We have applied this method to downloaded 
light curves with different binning (20 and 40 
days), and estimated ag i.e. intrinsic variabil- 
ity as well as variability strength Sv for all ASM 
detected AGN. This exercise is carried out for 
sum band light curve over 1.5-12 keV, A band 
i.e. 1.5-3 keV, as well as (B+C) band i.e. 3-12 
keV. Variability strengths for 20 and 40 days bin- 
ning were found to be similar. Results for 20 days 
binning for sum band data are given in Table [1] 
for 30 AGN which are common between ASM and 
BAT database. We have also listed ASM vari- 
ability strengths for four AGNs with ASM count 
rates above 0.2 c/s, but not detected by BAT. Also 
included in the table is BAT variability strength 
Mrk 3, which is not detected by ASM. In view 
of unknown systematics in ASM data, we have 
quoted ASM variability strengths in various bands 
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for only sources with count rates above 0.0895 
c/s for S band and 0.0274 c/s for A band. BAT 
variability strengt hs given in this table are from 
IBeckmann et al.l ( 20071) . We have quoted their 
results as upper limits for the sources where their 
estimates of variability strength are negative or 
smaller than error on estimate. 

Since the X-ray spectral parameters (like the 
power-law index) are not available to all these 
sources, we have calculated the X-ray luminosi- 
ties by converting the ASM sum band rates (R 
counts/s) to to energy flux using 

F[erg/crr? /s] = 3.2 x 10"^" x R[cts/s] (12) 

Th is assumes Cr a b-like spectrum, as prescribed 
bv iGrimm et al. ( 2002h . It should be noted that 
this assumption is not strictly valid for absorbed 
sources. For these sources photon index, over the 
A band in particular, could be much flatter than 
that of Crab. For the sources where spectral pa- 
rameters are available, we get an average conver- 
sion factor very close to this value (3.0 x 10~^° x 
R [cts/s]). 

To compare the flux in Crab units we have 
used a conversion factor of 75.5 counts/s (ASM 
sum ban d) and 453.8 coun t s/s for BAT rates as 



indicated b y much lower variability strength for 



this source (jSanders fc: Fabian 1 120071 ) . Solid hue 



given by IBeckmann et al." ( 2007t ). This allows 



easy comparison of fluxes in two different energy 
bands. Fig. [S] shows variation of BAT flux vs 
ASM total flux in the energy range of 1.5-12 keV. 
Blazars are indicated with plus signs, Seyfert 1 
with asterisks, Seyfert 1.5 with diamonds, Seyfert 
1.8 and 1.9 by square, Seyfert 2 with triangles and 
radio galaxy with a cross. Error bars are indicated 
for both BAT and ASM fluxes. Typically these er- 
ror bars are smaller than symbol size. One Blazar, 
Mrk 421 and one Seyfert 2 galaxy, NGC1275 (not 
shown in the figure) seem to have much higher 
flux in ASM compared to other AGN with sim- 
ilar flux values in BAT energy band. Mrk 421 
is highly variable source as can be seen from the 
ASM light curve. BAT data corresponds to only 9 
months out of 12.5 years of data accumulated by 
ASM. (Mrk 421 was in a somewhat low state dur- 
ing the BAT observations.) In case of NGC 1275, 
ASM flux is about 33.7 mCrab whereas BAT flux 
is 4.5 mCrab. The field of view of NGC 1275 con- 
tains the Perseus cluster and X-ray emission from 
this cluster could be responsible for higher count 
rate and dilution of observed X-ray variability as 



in the figure indicates both fluxes being linearly 
related (as will be the case if the spectral slope of 
AGNs is crab like - a photon index of 2). BAT flux 
seems to be increasing rapidly relative to the ASM 
flux. The lack of a strong correlation between 
ASM and BAT fluxes could be due to the mixing of 
different types of objects: sources with very soft 
spectrum which could be preferentially observed 
by the ROSAT satellite do not show flux corre- 
lation between the soft and hard bands (Tueller, 
Mushotzky, Barthelmy et al. 2008), whereas heav- 
ily obscured sources are obviously are not expected 
to show correlation (because the very soft flux in 
these sources may not be originating from the cen- 
tral engine of the AGN). 

Fig. do] show s luminosity of AGN fr om BAT 
data as given bv lBeckmann et al. ( 20071 ) vs ASM 
luminosity. On logarithmic scale BAT luminosity 
seems to be increasing with ASM luminosity lin- 
early. There are two AGN which are slightly away 
from the trend shown by other AGN. These are 
Blazar : Mrk 421 and Seyfert 2 : NGC 1275 as in 
the case of Fig. |S| 

Fig. [n] shows the plot of ASM and BAT vari- 
ability strength in different energy ranges. Top 
panel shows ASM variability strength in 1.2-3 
keV band, second panel shows ASM variability 
strength in 3-12 keV band and third panel shows 
BAT variability strength in 14-195 keV as a func- 
tion of ASM variability strength in 1.5-12 keV 
range. Bottom panel also shows BAT variability 
strength in 14-195 keV as a function of ASM vari- 
ability strength in 1.5-12 keV range, for near si- 
multaneous data. Here we have selected ASM data 
corresponding to first nine months of BAT oper- 
ation. Error bars on ASM variability strength in 
this panel are large and for some of the sources 
only upper limits on variability strength could be 
estimated for ASM data. For the sake of clarity 
we have restricted X-axis range to 120 for all the 
panels. Bottom-most panel has three points with 
ASM variability strength exceeding this value. 
Solid line in each panel corresponds to slope 1. 
It can be seen that variability strength decreases 
with increase in energy from 1.2-3 keV to 14-195 
keV. Variability strength in 3-12 keV band seems 
to be well correlated and almost similar to the one 
in 1.5-12 keV band indicating that in these two 
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ASM flux {1.5-12 keV} in mCrob 

Fig. 9.— BAT flux (14-195 keV) vs. ASM flux 
(1.5-12 keV) for a sample of 30 AGN common in 
BAT and ASM data base. Blazars are indicated 
with +, Seyfert 1 by asterisk, Seyfert 1.5 by di- 
amond, Seyfert 1.8 and 1.9 by square, Seyfert 2 
by triangle and radio galaxy by cross. Solid line 
indicates slope 1. 



+ Mkn 421 
A NGC 1275 



oy I x(l,5-12 keV) erg/s 

Fig. 10. — Hard X-ray luminosity of AGN from 
BAT vs soft X-ray luminosity from ASM for 30 
common AGN. Blazars are indicated with 
Seyfert 1 by asterisk, Seyfert 1.5 by diamond, 
Seyfert 1.8 and 1.9 by square, Seyfert 2 by triangle 
and radio galaxy by cross. 



energy bands most of the flux is in the large over- 
lapping energy range. To investigate this further 
we have grouped the sources in three groups, 1. 3 
Blazars and one radio galaxy, 2. 7 Seyfert 1 and 4 
Seyfert 1.5 galaxies and 3. 1 Seyfert 1.8, 2 Seyfert 
1.9 and 11 Seyfert 2 galaxies. Fig. [T2l shows dis- 
tributions of variability strengths in 1.5-12 keV, 
1.5-3 keV, 3-12 keV and 14-195 keV bands. Dot- 
ted line corresponds histogram generated includ- 
ing all sources, dashed line corresponds to blazars 
i.e group 1 and solid line corresponds to Seyfert 1 
and 1.5, i.e. group 2. This figure includes BAT 
variability strengths given as upper limits in Ta- 
ble [H Average values of variability strengths for 
different groups and for different energy ranges 
are listed in Table [2l Values given here are the 
weighted mean and estimated error on the mean. 
BAT variability strengths given as upper limits in 
Table [1] are not included here. NGC 1275 is ex- 
cluded here for ASM in these calculations since 
variability strength for this object quoted here 
could be grossly underestimated due to contam- 
ination caused by Perseus cluster in the field of 
view of ASM. This table clearly shows the trend 
of decrease in variability strength with increase in 
energy for all types of AGN. 

Higher variability strength in 1.5-3 keV band 
compared to the one in 3-12 keV band for 
Sey fert I's is cons istent with previous findings 
fsee lArevalo et al"~1 12008) and references therein). 
They have found that the variability as a function 
of energy peaks around 2 keV one time scales of 
one to a few days. Our result extends this prop- 
erty to longer time scales. 

We have compared our results with the vari- 
ability stren gth in 2-12 keV from RXTE/P CA as 
presented bv iMarkowitz fc Edelson I (|2004[ ) . They 
have given variability strength for several Seyfert 1 
and 1.5 galaxies on different time scales. There are 
five sources common between their sample with 
timescale of 1296 days (Table 3 in their paper) 
and the present sample. These are 3C120, NGC 
3783, NGC 3516, MCG-6-30-15 and NGC 3227. 
For these sources we have computed average vari- 
ability strength in the energy range of 1.5 - 12 keV 
with 40 days binning covering rough l y the same 
MJD range as IMarkowitz fc Edelsonl (|2004l ). Bin 
size of 40 days used here is close to 34.4 days bin 
size used by them. Fig [12] shows the comparison 
of ASM variability strengths in 1.5 - 12 keV with 
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Table 2 

Average variability strengths in different energy ranges 



Energy Range All sources Blazars Sy 1-1.5 Sy 1.8-2 



ASM (1.5-12 keV) 41.1±0.8 

ASM (1.5-3 keV) 61.4±1.5 

ASM (3-12 keV) 44.3±1.0 

BAT (14-195 keV) 12.6±1.3 



57.3±1.7 38.4±1.4 34.1±1.2 

64.7±2.0 64.0±2.9 51.5±3.0 

64.2±1.9 39.1±1.9 36.4±1.4 

17.2±3.6 19.2±5.3 11.4±1.5 




Fig. 11. — (a) ASM variability strength in 1.5- 
3 keV band, (b) ASM variability strength in 3- 
12 keV band (c) BAT variability strength in 14- 
195 kcV band as a function of ASM variability 
strength in 1.5-12 keV band and (d) BAT variabil- 
ity strength in 14-195 keV band as a function of 
ASM variability in 1.5-12 keV band for simultane- 
ous data. Variability strength decreases at higher 
energies. 



Fig. 12. — Distribution of variability strengths of 
AGN in (a) 1.5-12 keV, (b) 1.5-3 keV, (c) 3-12 
keV and (d) 14-195 keV as obtained from ASM 
and BAT data. In each panel, dotted line repre- 
sents distribution including all objects, dashed line 
corresponds to blazars and radio galaxy and solid 
line represents variability strength distribution for 
Seyfert 1 and 1.5. 
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PCA variability in 2 - 12 keV range. ASM vari- 
ability seems to be higher than the PCA variabil- 
ity. 

3.4. Variability strength vs break time 
scale 

AGN variability is often expressed in terms of 
fluctuations in power spectral density (PSD) i.e. 
variability power P{i') as a function of frequency, 
v. On longer timescales, the PSDs of AGN are 
fitted by a power-law of slope -1 which breaks to a 
steeper slope (> 2) at ti mescales shorter than th e 



break' timescale, Tg (jMarkowitz et al. I l2003l) 



Break timescale is expected to depend on the black 
hole mass. Considering similarities between PSDs 
of X-ray binaries and AGN, break timescales for 
AGN are expected to be in the range of < 1 day 
to > a hundred days. Earlier attempts of com- 
paring the Tb with black hole mass {Mbh) for 
AGN ha ve shown a rough linear s caling but with a 



scatter ( Markowitz et al. 2003t McHardv et al 



20041 ). Improving on it, an inverse depen- 
dence on a second variab le, probably the accre- 
tion rate, was suggested ( McHardv et al. 2004 



20051: Uttlev &: McHardv Il2005^ . iMcHardv et al . 
( 20061 ) quantified the relationship between Tb, 
yLBH and Li,oi (in place of the accretion rate), as 
logTB = AlogMBH - BhgLboi + C, where Mbh 
is in units of 10® solar masses, Lboi is bolometric 
luminosity in lO"*"* erg s~^. Best fit values deter- 
mined by them are A (= 2.1), B (= 0.98) and C 
(= 2.32). Using this relation, we have calculated 
break timescales (Tb) for our AGN sample. Values 
oflog(MBH) (in units of IO^Mq) and log(Lf,oi) (in 
units of lO"^"^ erg s~^ ) listed in Table [T] are adopted 
from I Woo fc Urrv I liool. lUttlev fc McHardv 



(l2005f ) and IWandel & Miishotzkv I (|l98ffl . The 

masses are estimated either using reverberation 
mapping or the BLR size-luminosity relation or 
the stellar velocity dispersion method. Lboi values 
are also taken from references mentioned above. 
It should be noted that Lboi for the blazars may 
be beamed, leading to a wrong estimates of Tb ■ 

Figure [14] shows ASM variability on the scale 
of 20 days as a function of log(rB). Figure [T51 
shows BAT variability strength as a function of 
log(TB). Out of 44 AG N from BAT sample of 
Beckmann et al. ] (|2007l ). black hole masses and 



be estimated for these sources. Out of these, 
19 sources are common between ASM and BAT 
database. The remaining source is Mrk 3, which 
is not in the ASM source list. Figure shows some 
increase in variability of BAT data for AGN in 
the neighborhood of 20 days and decrease on both 
lower and higher sides. The average variability for 
log(TB) between -1 to -1-1.6 is 62.9 and 53.6, for the 
ASM and BAT respectively and these values are 
27.4 and 16.7 outside this range. The decrease for 
higher Tb is as expected. The lower variability at 
lower Tb could be due to a variety of reasons like 
a) inclusion of a blazar for which the McHardy et 
al (2006) relation may not be valid b) at very low 
value of Tb, we are sampling the sources at very 
low frequencies where a second turn-over of PSD 
is likely if the sources are a state like the low-hard 
states of Galactic black hole sources. 

3.5. Energy spectra 

We checked for data from pointed observations 
of these AGN by X-Ray Telescope (XRT) onboard 
Swift during this 9 month period. Data was avail- 
able for four sources from our list. These were 
three Blazars : 3C273, 3C454.3 and Mrk 421 and 
one Seyfert 2 Galaxy : Gen A. We have fitted 
spectrum over the energy range of 0.3-10 keV for 
these four objects. For each of these objects source 
and background photons were selected using XS- 
ELECT version 2.4. For Photon Counting (PC) 
mode source photons were selected in a circular 
region with the radius of 20 pixels (i.e 47 arc- 
seconds), whereas background photons were ex- 
tracted from nearby circular region with a radius 
of 40 pixels. For data collected in Windowed Tim- 
ing (WT) mode, source photons were extracted 
using box region with the length of 40 pixels and 
width about 20 pixels. Events with grades 0-12 
and 0-2 were selected for PC and WT mode data, 
respectively. The spectral data were rebinned by 
GRPPHA 3.0.0 with 20 photons per bin. Standard 
auxiliary response files and response matrices were 
used. 

Spectra for these sources were fitted using 
XSPEC version 12.3.1 with a model consisting 
of absorbed power law over the energy range of 
0.3-10 keV. In each ca se Nh was fixed to th e 



value given in Table 2 of lBeckmann et al. I (|2007l ) . 



bolometric luminosity are available for 20 sources 
(jUttlev fc McHardv iBooih . and hence Tb could 



Power law indices and 2-10 keV flux obtained from 
fit are given in Tabled. It should be noted here. 
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Fig. 13.— ASM variability strength (1.5-12 keV) 
vs RXTE PCA variability strength (2-12 keV) for 
Seyfert 1 and Seyfert 1.5 galaxies. 
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Fig. 14. — ASM variability strength vs break time 
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Fig. 15. — BAT variability strength vs break time 
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however, that the Swift observations may not be 
equal to the average photon index representing 
the ASM observations. 

Among the three Blazars, Mrk 421 has the 
steepest spectrum and 3C 454.3 has the flattest 
one. In case of 3C454.3, we have not estimated 
the variability strength in ASM band because of 
low count rate. ASM variability strengths for 
other two blazars are : 83.7±4.2 for Mrk 421 
and 34.7±3.4 for 3C273. Whereas BAT variabil- 
ity strengths of these three sources arranged in the 
order of increasing power-law index are : 3C454.4 
(42±12), 3C273 (15±5) and Mrk 421 (142±38). 
This indicates an increase in variability strength 
for sources with steeper spectra. But this inference 
should be taken with caution. It should be noted 
that in case of Mrk 421, which is a BL Lac type 
of object, the hard X-ray band is dominated by 
high energy synchrotron emission whereas in case 
of 3C454.3 and 3C273, which are FSRQ, hard X- 
ray emission is mainly inverse-Compton emission. 
This could be the one possible reason for steeper 
spectrum of Mrk 421 compared to the other two 
objects. It is also quite possible that the long dura- 
tion variabilities of blazars are influenced by infre- 
quent strong flare events and the higher variability 
in any one blazar could simply be the effect of a 
few large events which happen to occur in that 
source. 

4. Discussion and conclusions 

Study of X-ray variability of AGN at longer 

time scales and also at different energies is very 
important to pin down the accretion disk geom- 
etry and the radiation processes involved in the 
X-ray emission very close to the black holes. We 
have made a systematic study of the soft X-ray 
variability characteristics of all AGN with mea- 
sured long time scale variability from Swift/BAT. 

One of the important findings of this work is 
that individual ASM dwells can be co-added to 
obtain flux integrated over long time scales and by 
propagating the measurement errors, very low er- 
rors on the data points are obtained. This method 
assumes that most of the systematic errors in the 
flux measurements in the individual dwell mea- 
surements are understood and taken care of. Sup- 
port for this assumption comes from the fact that 
the long duration light curves of several AGN bear 
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Table 3 
Spectral Fits 



Source 


Observation 


Nh 


Powerlaw 


2-10 keV 


name 


duration 


1022cni-2 


index 


flux 




s 






erg/cm^/s 


3C 273 


2962.6 


0.0316 


1.79±0.030 


5.87x10-" 


3C 454.3 


13670.0 


0.0631 


1.32±0.02 


5.43x10-" 


Mrk 421 


6958.0 


0.001 


2.38±0.01 


8.34x10-" 


Con A 


7699.5 


12.59 


1.64±0.06 


1.99x10-1° 



a striking similarity to the light curves obtained 
from pointed RXTE/PCA observations. Some 
more work, however, need to be done to under- 
stand possible time dependent systematic errors 
so that a complete data set can be used to derive 
meaningful power spectral densities. 

Another important finding is the quite uniform 
variation of variability strength with energy for di- 
verse classes of objects, except for the blazar Mrk 
421. For this particular source, there is a marked 
increase in the variability as a function of energy. 
Such behavior has been ob served at shorter dura- 
tions. Horan et al. 1 (l2009l) made a detailed multi- 
wavelength study for a period of about 300 days 
and find the variability to increase across the full 
electro-magnetic band. In particular, they noticed 
a sharp increase of variability from soft X-rays to 
hard X-rays: the values for variability (F^q^) were 
26.9, 44.3, 52.9 and 99.3, respectively for Swift 
XRT, RXTE-ASM, RXTE-PCA and Swift-BAT 
data (for observation bins of a few tens of thou- 
sand seconds and duration in the range of 24 to 
256 days). The values that we obtained are 79 
(RXTE-ASM 1.5-3 keV), 88 (RXTE-ASM 3 - 
12 keV) for 11 years and 142 (Swift-BAT for 300 
days). The energy dependent variability behavior 
of the the blazar in our list, 3C 273, for which data 
is available both in the soft and hard X-rays, shows 
an energy dependent behavi or similar to Seyfert 
1 galaxies. McHardv ( 2006f ) has noted that the 
wide band PSD of 3C273 is identical to that of 
Seyfert galaxies and contended that the process re- 
sponsible for variability in this source is the same 
that produces variability in non-beamed sources. 

Most of the other sources in our sample are 
Seyfert galaxies and a decrease in the variability 



as a function of energy is observed. In a detailed 
study of a large number of Seyfert galaxies based 
on the first seven years of RXTE -PCA monitor- 
ing, Markowitz fc Edelson ( 20041 ) give energy de- 
pendent variability data for long durations and a 
similar trend of decrease in variability with energy 
has been noticed. There are five Seyfert 1 galaxies 
(3C120, 3C390.3, IC4329A, NGC3783, and MCG- 
6-30-15) and two Seyfert 1.5 galaxies (NGC4151 
and NGC3227) common between their sample and 
our work. The ratio of variability for the soft (2 
- 4 keV) and the hard (7 - 12 keV) band are, 
respectively, 1.15, 1.26, 1.19, 1.28, 1.37, 1.26 and 
1.35, for a duration of 216 days. For compari- 
son similar ratios from our data (1.5 - 3 keV vs 
3-12 keV), for the above sources, respectively, 
are, 1.65, 1.62, 2.65, 2.9, 1.71, 2.6, and 2.5, for 
the full ASM duration. The typical uncertainties 
are 0.06 for the RXTE data and 0.3 for our data. 
Though this ratio from RXTE data is fairly same 
for all the sources (about 1.27), for the ASM data 
it varies from 1.66 to 2.64. For five sources en- 
ergy dependent variabili ty at longer times scales 
(> 16 days) is given bv iMarkowitz et al. I (|2003l ) 
and for two of the sources where we have ob- 
tained very high ratio (NGC 3783 and NGC 4151) 
a sharp upturn towards lower energies are noticed 
(see their Fig 6). Further, RXTE-ASM is more 
sensitive to low energy X-rays. The effective area 
at 2 keV is about 33 % of the effective area at 5 keV 
( Levine et al. Ill996l). whereas this nu mber is < 1% 
for RXTE PGA (jjahoda et al. llioo e'). mainly due 
to the propane layer on the top side of RXTE-PGA 
which has lower transparency at low energies due 
to aluminized mylar windows and the permeation 
of xenon gas into the propane layer (|Jahoda et al~l 
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19961) . Hence it is quite possible that the higher 
variabihty at low energies derived by us is real and 
indicate a sharp increase of variability in AGNs at 
lower energies for longer time scales. Though we 
cannot completely rule out some further system- 
atic errors, it is quite evident that RXTE-ASM 
is an unique instrument to probe long time scale 
variability at lower energies. 

There are, indeed, several indicators from early 
EXO SAT data which supports th e above asser- 
tion. IWalter fc C ourvoisier I ( 1992f ) discuss X-ray 
variability of several AGNs observed by EXOSAT 
over the period of an year and note that soft X- 
rays vary by factors up to 7, much higher than that 
seen in hard X-rays. Two of the sources in their 
list is common to the present work: MR2251-178 
and NGC 3783. For these two sources, they report 
a peak-to-peak variation of 4.5 and 3, respectively 
in the soft X-rays, compared to 1.9 and 1.8 in hard 
X-rays. This should be compared with the Sy de- 
rived by us for these two sources: 97 and 67, re- 
spectively in the 1.5 - 3 keV band and 42±14 and 
24±7 in the 3-12 keV band, respectively. Though 
there is a strong indication of very high soft X-ray 
variability at long time scales, a detailed compar- 
ison with a large number of pointed observations 
for individual sources are required to pin down 
any residual systematic errors in the RXTE-ASM 
data. 

Energy dependent variability has been stud- 
ied extens i vely in Galactic black hole sources (see 
Zdziarski I (j2005l ) and references therein) and one 
of the models used to explain the energy depen- 
dency is the variations in the input parameters of 
the thermal Comptonization process which results 
in an pivoting power-law spectrum. Though the 
very low value of variability in the BAT energy 
band may be influenced by a non-varying spec- 
tral component in this energy band (like reflec- 
tion), detailed wide band spectroscopy of Seyfert 
galaxies suggest a common continuous phenomena 
for the decreasing variability with increasing en- 
ergy. For example, the radio galaxy 3C120 shows 
only a moderate reflection component and very 
strong spectral variability (over two days) with 
the spectral index correlated with the soft X-ray 
flux, indicating a phenomena where the input seed 
photon variation for a thermal Comptonization 
process causes the energy depende nt X-ray vari- 
ability (jZdziarski fc Grandi M200lh . If the same 



phenomena is responsible for most of the sources 
reported in this work, it suggests that a) for a ma- 
jority of AGNs the thermal Comptonization is the 
dominant process, b) the pivot energy is higher 
(greater than about 20 keV). 

This work is based on results provided by the 
ASM/RXTE teams at MIT and at the RXTE SOF 
and GOF at NASA's GSFC. This research has 
also made use of data obtained from the High En- 
ergy Astrophysics Science Archive Research Cen- 
ter (HEASARC), provided by NASA's Goddard 
Space Flight Center. We are very grateful to the 
referee of this paper for the very thoughtful and 
critical comments. 
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